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000000 TITTTITITT SSSSSSSS_ PPPPPPPP 000000 WW WwW RRRRRRRR UU UU 
000000 TITTTTTITT SSSSSSSS_ PPPPPPPP 000000 wis WwW RRRRRRRR 

00 00 TT SS PP 00 00 WwW 

00 00 TT $$ PP PP 00 00 WW Ww RR RR UU UU 

00 00 TT $$ PP PP 00 00 WW Ww RR RR UU UU 

00 00 TT SS PP PP 00 00 WwW WW RR RR LL UU UU 

00 00 TT SSSSSS PPPPPPPP 00 00 WW Ww RRRRRRRR LL UU UU 

00 00 TT SSSSSS PPPPPPPP 00 00 WW WW RRRRRRRR LL UU UU 

00 00 TT SS PP 00 00 WW WW WW RR RR LL UU UU 

00 00 TT SS PP 00 00 WW WW Ww RR RR LL UU UU 

00 00 TT SS PP 00 00 WWwWW WwWWw RR RR LL UU UU eeee 

00 00 TT SS PP 00 00 WWwWW WWWWw RR RR LL UU U eoee 
000000 TT SSSSSSSS PP 000000 wil Ww RR RR LLELLLLLLLL 9 UUUUUUUUUU cece 
000000 TT SSSSSSSS PP 000000 wid Ww RR RR LLLLLLLLLL 9 UUUUUUUUUU cece 

LL III! SSSSSSSS 

LL HII] SSSSSSSS 

LL I] SS 

LL I] SS 

LL I] SS 

LL I] $$ 

LL I] SSSSSS 

LL I] SSSSSS 

LL I] SS 

LL I] SS 

LL I] SS 

LL I] SS 

LLELLLLLLLL III] SSSSSSSS 

LLLLLLLLLL HII] SSSSSSSS 
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OTSSPOWRLU = F_floating ** unsigned integer power routine 
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COPYRIGHT (c) 1978, 1980, 1982, 1984 B 
DIGITAL EQUIPMENT CORPORATION, MAYNARD. MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


*® ® 
® ® 
*® ® 
® ® 
:@ ® 
: OFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED * 
ie NoyAGCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
:e ION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER 
ie THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY © 
ie PER O TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 
:® TRANSFERRED. . 
:@ * 
;® THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE * 
® AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT * 
:® CORPORATION. * 
:* * 
‘* ® 
‘* * 
:¢ a 
*® * 
® ® 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 
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-+¢ 
; FACILITY: Language support procedures, Mathematics division 
; ABSTRACT: 


yet module contains OTSSPOWRLU, a procedure which raises an 
F_floating value to an unsigned integer power. 


ENVIRONMENT: Runs at any access mode, AST Reentrant 
; AUTHOR: Steven B. Lionel, CREATION DATE: 22-JUL-1981 
; MODIFIED BY: 
| 1-001 = Original. SBL 22-JUL-1981 | 
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7 ~SBTTL DECLARATIONS 
4 : LIBRARY MACRO CALLS: 
0° SSSDEF : System error codes 
1 SCHF DEF 3; Define condition handler symbols. 
§ SSFDEF ; Define stack frame symbols. 
SPSLDEF ; Define program status longword 
¢ . 3; symbols. 
$ ; EXTERNAL DECLARATIONS: | 
8 : -DSABL GBL ; Force all external symbols to be declared 
00 -EXTRN MTHSK_UNDEXP ; Undefined exponentiation 
$8 60 -EXTRN MTHSK_FLOOVEMAT ; Floating overflow in Math Library 
0 61 -EXTRN MTHSK_FLOUNDMAT ; Floating underflow in Math Library 
9 66 .EXTRN MTHSSSIGNAL ; Math error routine 
0000 64 ; MACROS: | 
000 65 ; 
44 $6 3 NONE 
000 67 ; 
000 68 ; EQUATED SYMBOLS: 
000 9; 
000 0; NONE 
444 71 ; 
00 72 ; OWN STORAGE: 
000 Ts 2 
000 74; NONE 
000 2 3 
44 16 ; PSECT DECLARATIONS: 
00000009 78° -PSECT _OTSSCODE PIC, USR, CON, REL, LCL, SHR, - 
9000 43 EXE, RD, NOWRT, LONG 
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M88 OTSSPOURLU'S F floating ae cattonod inte S-sePoloRe STiobsco EMTMNTE eee vorsecumcumaR:?  o% (dy 
is -SBTTL OTSSPOWRLU = F_floating ** unsigned integer power routine 
; FUNCTIONAL DESCRIPTION: 


Raise an F bag thee: value to an unsigned integer power giving 
an F_floating result. 


The result is given by: 


J 
m” 
4 ; 
: 
0 : : 
p 31 : base exponent result | 
38 : any >0 product (base * 2**i) where i is each 
Be : non-zero bit position in exponent 
0 96 : > = 0 1.0 | 
0 97 ; = 2 Undefined exponentation 
000 98 ; < = 0 1.0 
000 99 ; 
44 100 ; 
000 101 ; CALLING SEQUENCE: 
000 196 : 
4 : Z : result.wf.v = OTSSPOWRLU (base.rf.v. exponent.rlu.v) 
000 198 : FORMAL PARAMETERS: 
000 196 ; 
000 10 
00000004 44 13 base = 4 ; F_floating base 
00000008 44 419 exponent = 8 ; Unsigned longword integer exponent 
000 We ; 
000 1135 ; IMPLICIT INPUTS: 
000 114; 
4 115 ; NONE 
00 16 3 
000 117 ; IMPLICIT OUTPUTS: 
0000 118 ; 
000 119 ; NONE 
00 120 ; 
bo : 1 ; ROUTINE VALUE: 
bp : : The base raised to the exponent's power. 
99 5 ; SIDE EFFECTS: 
0 1 § : May signal: 
0 1 3 3 MTHS_FLOOVEMAT - floating overflow in Math Library 
0 129 ; MTH$S_FLOUNDMAT - floating underflow in Math Library, only 
0 130 ; if the caller has FU enabled. 
ta) 3 MTHS_UNDEXP - undefined exponentiation, if both base and 
: § : exponent are zero. 
4 to- | 
0004 p : § ENTRY OTSSPOWRLU, “M<R2> ; Entry point 
64°AF = 9E 138 MOVAB B*EXC_HANDLER, (FP) ; Translate exceptions to MTH$ errors 
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40 oF 6 139 BISPSW #PSLSM_FU ; Enable floating underflow detection 
8 A 179 MOVF wi, R ; RO = initial result 
2] 4 AC D «14 MOVF base(AP), R1 ; R1 = base 
2 08 at D 146 MOVL exponent (AP) , R2 : R2 = exponent 
1 1 14 BNEQ BEGIN : frie if not zero 
oS 1 144 TSTF & ; Is base also zero? 
pa 122 BEQL UNDEF INED ; jf 30. error. Otherwise, answer is 
04 0018 123 RET : with result of 1.0 
1€ 409 ss 
1¢ 139 : Do the first iteration here so that we can clear the high-order exponent 
Bie 1 ; Ha Afterwards, that bit will be zero so we don't have to worry about 
ie aa 
01C 136 BEGIN 
03 5 00 4 01 138 BBCC #0, R2, 10% ; Skip initial multiply 
5 51 0 00 15 MOVF R1 ; New partial result 
 —_—-. fa - sem 83 138 10$: ROTL #-1, R2, R2 3 eye, hey tow exponent bit and clear 
10 13 00 ; 160 BEQL DONE : Is-that all? If so, we're done. 
ee 
$0 A 188 : Each time we get here we know that there's at least one more exponent bit 
O02A 164; left, so square our current power of the base. 
Q02A 165 ;:- 
ae 
31hUC«<C KCC 8 A 168 MULF2 R11, R1 ; Square base. This may overflow or 
002D 169 ; underflow. In either case, the 
88 D 170 : final result would also overflow 
02D 171 : or underflow. Our exception 
Bp D 78 ; handler will catch these cases. 
03 52 €9 O002D 17 BLBC R2, NEXT_BIT ; See if this exponent bit is on. 
BS 0 174 : If so, multiply base by product. 
50. 651 44 00 175 MULF2 R1, RO ; This too may overflow or underflow, 
0 176 3 which will be caught by the handler. 
0 177 NEXT_BIT: 
52 52 FF BF 8 0 178 ASHL #-1, R2, R2 ; Get next bit of exponent. 
FO. 0 : 179 BNEQ LOOP : Keep going if not all done. 
O3A 181 DONE: 
04 8 ean ¢ RET ; Return with answer in RO. 
038 i 4 3+ 
6 185 : Undefined exponentation error - 0**0 
189 
0 1 8 UNDEFINED: : : : 
OO'8F 9A 1 MOVZBL #MTHSK va * wh -(SP) ; Indicate undefined exponentiation. 
50 O01 OF 78 F 190 ASHL , aT, R ; RO = reserved floating operand 
15 43. (191 BRB SIGNAL_ERROR 
28 198: 
43 194 : EXC_HANDLER unwinds here if we got a floating overflow. 
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OVERFLOW: 
7E  00°8F 3A MOVZBL #MTHSK a ote: -(SP) ; Floating overflow in Math Library 
50 «(01 F 8 ASHL #15, #T, R ; RO = reserved floating operand 
if BRB SIGNAL_ERROR 


+ 
EXC_HANDLER unwinds here if we got a floating underflow. 
If our caller enabled FU, then signal FLOUNDMAT, with a result of zero. 
Otherwise quietly return zero. 


ce 


NDERFLOW: 
CLRF RO i Initial result 
BBC #PSLSV_FU, SFSW_SAVE_PSW(FP), DONE ; Just return if 


6 Cc Ww. 
MOVZBL #MTHSK_FLOUNDMAT, -(SP) ; Floating underflow in Math Library 
BRB SIGNAL _ERROR 


50 
E404 AD 06 €1 
7E O0°SF 9A 


+ 


Signal a MTH$ error whose code has already been pushed on the stack. 
RO has either a reserved operand or a zero, depending on the exception type. 


GP te ee eee 


IGNAL ERROR: 
LRL ; Cancel our condition handler 
6D CLR (FP) C l diti handl 
Q00000000'GF 01 FB CALLS #1, G*MTHSSSIGNAL 3 convert to 32-bit condition code 
; and signal a MTH$ error. Omit 
3; second argument to show that this 
; is a CALL entry point 


RET ; Return to caller with result in RO 
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4 9 ;+ 
e4 ; Exe HANDLER = This condition handler gets control if we got an exception 
64 ; while performing the exponentiation. If it is a floaging overflow, 
64 ; then continue execution at OVERFLOW, to cause MTH$_FLOOVEMAT to be 
64 3 pignat ied. Similarly, continue at UNDERFLOW if an underflow was seen. 
64 4 ; Otherwise resignal. 
64 5 ;- 
64 § 
64 EXC_HANDLER: 
0000 2 8 «WORD “M<> ; Entry point 
50 4 + D 28 rt MOVL CHFSL_MCHARGLST(AP), RO ; Get mechanism arguments List 
A D A 4 TSTL CHFSL_MCH_DEPTH(RO) ; At depth zero? 
244—Ci«d1 6D 4g BNEQ RESIGNAL ; If not, resigna 
50 O04 AC OD oot 4 MOVL CHFSL_SIGARGLST(AP), RO ; Get signal arguments List 
51 04 A0 00 8 ; 44 MOVL CHFSL_SIG_NAME(RO), R1 ; Get signal name 
0484 8F 51 B1 7 45 CMPW R1, #SS$_FLTOVF_F : Overflow fault? 
18 13 007¢ 46 BEQaL «GoTo OVERFLOW ; If so, continue at OVERFLOW 
048¢ 8F 251~—soaB ore 4 CMPW ORT, BSS$_FLTOVE : Overflow trap? 
14 13 008 48 BEQL  GOTO_OVERFLOW : If so, continue at OVERFLOW 
04c4 8F 51 81 0085 49 CMPW = ORT, #SS$_FLTUND_F : Underflow fault? 
18 13 OO8A 50 BEQL GOTO_UNDERFLOW :; If so, continue at UNDERFLOW 
049C BF 651 B81 008C 51 CMPW R1, #SS$_FLTUND ; Underflow trap? 
11. 13 0091 3g BEQL § GOTO_UNDERFLOW ; If so, continue at UNDERFLOW 
ad 53 RESIGNAL: 4 
50 600918 8F =63C «(009 54 MOVZWL #SS$_RESIGNAL, RO ; Resignal exception : 
04 Bn88 3? ET ; Return to VMS exception dispatcher 
0099 3§ GOTO_OVERF LOW: 
51 60 00 0099 58 MOVL CHF$L_SIG_ARGS(RO), R1 ; Get number of signal arguments 
FC A041 AG AF DE 009C 259 MOVAL B*OVERFLOO, -4(RO)ERI] : Move address of OVERFLOW routine 
O0A 60 3; to PC in signal argument List 
09 11 Baa 9] BRB CONTINUE 3; Continue execution 
Q0A4 § GOTO_UNDERF LOW: 
51 60 DO Q0A4 64 MOVL CHFS$L_SIG_ARGS(RO), R1_ ; Get number of signal arguments. 
FC A041 =60AS AF DE OOA7 = 265 MOVAL B*UNDERFLOW, -4(ROSCRI] ; Move address of UNDERFLOW routine 
OAD 96 3 to PC in signal argument List 
ae rt 3 BRB CONTINUE ; Continue execution 
OAD 269 CONTINUE: 
50 01 3¢ OOAD 70 MOVZWL #SS$_CONTINUE, RO ; Continue at PC in signal List 
04 0981 re RET ; Return to VMS exception dispatcher 
1 7 
0081 74 - END : End of module OTS$POWRLU 
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Allocation PSECT No. Attributes 


00000000 ( 00 ba NOPIC USR CON ABS LCL 
0000000 ¢ 3 01 3°} NOPIC USR CON ABS LCL 
0000081 ( 17 02 2.) PIC USR CON REL LCL 


weer een en soso scene een eee a } 


“rts NOEXE NORD NOWRT NOVEC BYTE 
SHR EXE RD WRT NOVEC BYTE 
SHR EXE RD NOWRT NOVEC LONG 
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Page faults CPU Time Elapsed Time 
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ihe yorking (83 Limit was 1050 pages. 
3 ~ + ( popes) of virtual memory were used to buffer the intermediate code. 

seal were 30 pages of symbol table space allocated to hold 508 non-local i local symbols. 

74 source tinge were read in Pass 1, producing 13 object records in Pass 2. 

1 pages of virtual memory were used to define 10 macros. 
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! Macro Library statistics ! 
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Macro Library name Macros defined 
"$25SSDUAZ8:(SYSLIBISTARLET.MLB;2 (<ti‘(sté‘<t~S;*”S”S™S . 
562 GETS were required to define 7 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:OTSPOWRLU/OBJ=OBJ$:OTSPOWRLU MSRC$:OTSPOWRLU/UPDATE=(ENH$:OTSPOWRLU) 
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